NASA Contractor Report 4612 



On the Various Forms of the Energy Equation 
for a Dilute, Monatomic Mixture of Nonreacting 
Gases 

Christopher A. Kennedy 

7-/7 


(NASA-CR-4612) ON THE VARIOUS 

forms of the energy equation for a 

DILUTE, MONATOMIC MIXTURE OF 
NONREACTING GASES Final Report 
(Cal i forni a Uni v. ) 19 p 


N94-3 5864 
Unci as 

Hl/02 0016064 


Grant NAG1-1193 
Prepared for Langley Research Center 


July 1994 



NASA Contractor Report 4612 



On the Various Forms of the Energy Equation 
for a Dilute, Monatomic Mixture of Nonreacting 
Gases 


Christopher A. Kennedy 

University of California at San Diego • La folia, California 


National Aeronautics and Space Administration 
Langley Research Center • Hampton, Virginia 23681-0001 


Prepared for Langley Research Center 
under Grant NAG1-1193 


July 1994 




ON THE VARIOUS FORMS OF THE ENERGY EQUATION 
FOR A DILUTE, MONATOMIC MIXTURE 
OF NONREACTING GASES 


Christopher A. Kennedy 

Department of Applied Mechanics and Engineering Sciences 
University of California, San Diego, La Jolla, California 92093 


Introduction 

Choosing the appropriate form of the governing equation set to solve in 
fluid dynamics is the first step in the solution of a particular problem. In 
the case of gas mixtures, the governing equations become rather formidable 
and a complete listing of the equations in their various forms and methods 
to evaluate the transport coefficients is difficult to find. This paper seeks to 
compile common as well as less well known results in a single document. 

Forms for the various quantities involved in multicomponent hydrodynam- 
ics may be derived either phenomenologically 6 ’ 18, 10 or by methods of kinetic 
theory 31 ’ 11 , Ideally, the full limitations of the derivation need to be stated a 
priori. For this reason, we choose the kinetic theory approach. In addition, 
kinetic theory allows us to explicitly calculate the relevant transport coeffi- 
cients. This paper will focus on the issue of complications to the governing 
equations because of multiple species but will not concern itself with issue of in- 
ternal degrees of freedom of the molecules leading to relaxation phenomena 34 
or antisymmetric stress 27, 24 , exothermic heat release that may cause high- 
energy “tails” to the distribution function of particle velocities 20 and differing 
temperatures of certain species 23 , or the effects of ternary, quartenary, etc. 
collisions found in dense gases 14 and liquids. Having made those severe restric- 
tions, we also limit discussion to small perturbations from equilibrium where 
one may safely assume that the first nonequilibrium expression derived from 
the Chapman-Enskog-Burnett 3 method of approximating the solution to the 
Boltzmann equation is valid. In terms of classical irreversible thermodynamics, 
“local” equilibrium prevails. This means that although the flow is not in equi- 
librium, locally, equilibrium thermodynamic variables are still defined. Results 
using the Chapman-Enskog-Burnett method are not valid for severe spatial 
gradients and for sufficiently high frequency phenomena. This amounts to no 
variation in quantities whose characteristic lengths and times are of the order 
of the mean free path or the time between collisions. Collision times are also 
presumed to be much smaller than the time between collisions. In cases where 
more severe phenomena occur, local equilibrium does not exist and unique def- 
initions for the nonequilibrium temperature and entropy are not forthcoming 2 . 
This is the topic of “Extended Irreversible Thermodynamics 26, 12 ”. Techniques 
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such as Grad’s moment method 9 , which allow for terms such as the heat flux 
vector, stress tensor, and species flux to become independent variables, are 
likely to lead to more accurate approximations to the Boltzmann equation in 
situations far from equilibrium 29 . Although the governing equations derived 
by Grad’s method are more complicated than those traditionally confronted, 
they are likely to be essential in describing the finer points of polyatomic and 
reacting dilute gases 33 as well as rarefied gases 35 where the Knudsen number 
may no longer be considered a small expansion parameter. 

It is the goal of this paper to summarize the various relationships between 
equations describing conservation of energy for a dilute, monatomic, nonreact- 
ing gas in local equilibrium. The gas is treated as nonrelativistic, not subject 
to magnetic or electric fields, or radiative effects. 

Energy Equations 

Using the first nonequilibrium approximation according to the Chapman- 
Enskog-Burnett theory in the kinetic theory of gases 31, 8 (the Navier-Stokes 
level), 
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where i and j are species indicies, a and /? are spatial indicies, ncs is the 
number of chemical species, q a is the heat flux vector, q a ( ie a) is the reduced 
heat flux vector, is the diffusion velocity, <jp a is the stress tensor, of tQ is the 
diffusion driving force vector, Df is the thermal diffusion coefficient, Dff is the 
multicomponent diffusion coefficient, Ao is the partial thermal conductivity, A M 
is the second coefficient of viscosity, p is the shear viscosity, X{ is the mole 
fraction, Y, is the mass fraction, hi is the partial enthalpy per unit mass, p 
is the pressure, T is the temperature, p is the density, / 1Q is the body force 
per unit mass, 9 is the dilatation, S a p is the kronicker delta, and D a ^ is the 
rate of deformation tensor. Several comments about these terms are in order. 
There are two different definitions in the literature for the multicomponent 
diffusion and thermal diffusion coefficients, those of Waldmann 30 and those 
of Curtiss et. al. 4, 16 . It is essential that readers be aware of which one they 
are confronting at any given moment. We chose those of Waldmann for their 
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simplicity 22, 5 and their straightforward correspondence to coefficients found in 
linear irreversible thermodynamics 11, 15 . For the second coefficient of viscosity, 
Stokes hypothesis holds rigorously and we may use pg = (A„ + | p) = 0, where 
Pb is the bulk viscosity. Note that the partial thermal conductivity is not the 
conductivity that is measured in experiments because thermal diffusion effects 
can not be completely isolated. We also note that all diffusion velocities and 
diffusion driving force vectors are not independent but are related through 
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By defining the velocity gradient tensor as 

L a (S = V/3U a (7) 

the rate of deformation tensor and its deviator are given by 

1 0 1 
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where u a is the barycentric, hydrodynamic, or center-of-mass velocity of the 
gas. It will also be useful to write out the equations of state (9), species 
continuity (10), overall continuity (11), and momentum conservation (12) for 
future reference; 
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where the reaction rate, th;, is presumed to vanish and n and k are the number 
of atoms per unit volume and Boltzmann’s constant. The conservative forms 
of equations (10) and (12) may be obtained by simply adding zero to both 
equations in the form of Eq. (11) times Y{ and respectively. Use of the 
following relations has been made in order to transform the equation of state; 
n — n t -, pi = n t *m M p = Pi > where m t being the mass per atom, 

Yi = P J, Xi = R = % , R° = kN a , Wi = Nam,, where W, R°. N a are 
the atomic weight of species i, the average atomic weight of the mixture, the 
universal gas constant, and Avogadro’s number. Clearly Xi = ^"=1 Y, = 
1. 

The mole and mass fractions may be related by, 
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and the average molecular weight, W, is given by 

ncs 

= -£,(x t w,) 

k = 1 

These may be related through 

W L= Yi_ 

W X , 

Equations for the heat flux vector and diffusion vector appear in many 
different, but equivalent forms. Beginning with the expression 25, 17, 16 
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where Df- is the binary diffusion coefficient, we may rewrite the diffusion 
velocity equation (3) as the Stefan-Maxwell equation (17) using ~ 6 

and $Zr=i dia = 0 as 
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where the thermal diffusion ratios, kf, and antisymmetrical thermal diffusion 
factors, qT. , are defined by 
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Because there is no explicit mole fraction dependence in the thermal diffusion 
factors, they are sometimes the preferred term to describe thermal diffusion 
effects. There are [ncs(ncs + l)/2] independent multicomponent diffusion co- 
efficients and thermal (ncs — 1) independent thermal diffusion coefficients, as 
can be seen in the relations 
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Each of the terms, , D fj, and Rij = are symmetrical, the last 

term being referred to as the impedance matrix. In the special case of a 
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binary mixture, the multicomponent diffusion coefficients are related to the 
binary diffusion coefficient as 
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Solving the Stefan-Maxwell equation (17) for and then substituting the 
result into the original expression for the heat flux vector (1) gives 
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An experimenter measuring thermal conductivity is likely to measure it 
when all diffusion velocities vanish at steady state. For this reason the thermal 
conductivity is defined in terms of the partial thermal conductivity as 
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1 :=1 

so that in the absence of a diffusion velocity, the coefficient in front of X a T is 
the value which is measured. With this we may write four equivalent forms of 
the heat flux vector (1), 
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where | kT is the enthalpy per atom, or with the use of the equation of state, 
|£, and hi is the partial enthalpy per unit mass. The relation between the 
thermal and partial thermal conductivities may also be written in terms of the 
thermal diffusion coefficients instead of the thermal diffusion ratios as 
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The energy equation may be found in the literature in many forms. It may 
derived as 6, 18, 10 
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TOTAL ENERGY EQUATION (e 0 ) 
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Dotting the momentum equation (12) with the center-of-mass velocity, u a , 
yields 

MECHANICAL ENERGY EQUATION (|u a u a ) 
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With the thermodynamic relation h = e + one may obtain 
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By considering h = h(Yk,T,p) and e = e(Yk,T,p ) in flows where local 
equilibrium exists, we may write the total differentials of the enthalpy and 
internal energy as 32 
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These may be simplified using (§f) y = C p , the isobaric heat capacity, 

(^)y p = Cv, the isochoric heat capacity, (Jpj-) T = hk, the partial enthalpy 

per unit mass, (^')_, = hk — , where k is coefficient of isothermal 

compressibility and a is coefficient of thermal expansion at constant pressure, 
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perfect gas, a = ^ and k = ^ reducing equations (43) and (44) to 
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From equations (45), (46) and (35) we may write 
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which enables use to derive 
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With equations (36) and (47), one form of the temperature equation appears 
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By eliminating the material derivative of the pressure from equation (51) 
with equation (50) we have 

TEMPERATURE EQUATION 
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The expression for the time rate of change of entropy per unit mass may 
be written as 
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where J* is the flux of entropy and is written in the context of linear irreversible 
thermodynamics 10, 6 ’ 18 , with the use of p,- = A, — Tsi , as 
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where /q is the chemical potential per unit mass and s, is the partial entropy 
per unit mass. Therefore, 
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Equations (58) and (59) represent the reversible change of entropy. The ir- 
reversible component of the change in entropy is derived by noting that the 
entropy generation is of the bilinear form 

tr’ = ZJX 

= y, LXX ( 60 ) 


where J — LX , X , and L are the respective thermodynamic fluxes, thermo- 
dynamic forces, and Onsager coefficients. Summation is considered over all 
forces and fluxes of the same tensorial order. The four thermodynamic forces 
are the diffusion force (X-^), the thermal force (X^), and the scalar (X^ 0 )) 
and tensor viscous forces (X^) given by 7 
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and hence 
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The irreversible generation of entropy, <r s , is then given by 
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Because the second law of thermodynamics requires that entropy must increase 
until it is maximized at equilibrium, it is important to verify that cr a > 0. By 
manipulating the Stefan-Maxwell equation (17), using (20), and the definition 
of the impedance matrix, we may write the diffusion driving force vector as 
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and taking the reduced form of the heat flux vector (28), we have 
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Upon substitution of (79) and (80) into (78), all occurances of the thermal 
diffusion coefficients cancel, leaving 
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where the vanishing bulk viscosity has been dropped, and may be rewritten as 
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It is now clear that entropy production is never negative because A, \i and Df- > 
0 and all other terms occur in quadratic form (thermal diffusion coefficients 
may be positive or negative). 

Putting the irreversible entropy generation (77) and the entropy flux (58) 
into (55) yields 
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Each of the four fundamental transport coefficients appearing the equations 
(p, Ao, DfJ, Df) may calculated to any desired accuracy, and in particular 
those occuring in (16). Formally, the transport coefficients are calculated by 
evaluating the following matrix expressions 
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Formal expressions for the ( ncs x ncs) matricies H pq (p, q < 2) and L pq 
(p,q < 3) are given in the literature 1,19,13 . The variable £ denotes the or- 
der of the approximation and the vectors X, 8{, and A; are given by X = 

{ Ai, X2, * * • , X ncs } , 8 { = {^tl? ‘ • * > ncs } j A ,■ = {(6,1 — Tl), (^2 — Fl), • * # ncs“ 

F nc5 )}. The superscript J 1 denotes transpose. Note that the first approxima- 
tion to both the partial thermal conductivity and thermal diffusion coefficients 
vanish identically. 
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